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We consider the statics and dynamics of a flexible polymer confined between parallel plates both in the 
presence and absence of hydrodynamic interactions. The hydrodynamic interactions are described at the level 
of the fluctuating, compressible Navier-Stokcs equation. We consider two cases: (i) confinement for both the 
solvent and the polymer, and (ii) confinement for the polymer only (in a 3D solvent), which is experimentally 
feasible, for instance, by (optical) trapping. We find a continuous transition from 2D to 3D dynamic scaling as 
a function of decreasing degree of confinement within the de Gennes and the weak-confinement regimes. We 
demonstrate that, in the presence of hydrodynamics, the polymer's center-of-mass diffusion coefficient in the 
direction parallel to the walls scales differently as a function of the level of confinement in cases (i) and (ii). 
We also find that in the commonly used Langevin dynamics description, the polymer swells more parallel to 
the walls than in the presence of hydrodynamics, and the planar diffusion coefficient shows scaling behavior 
similar to case (ii) rather than case (i). In addition, we quantify the differences in the static structure factor 
of the polymer between cases (i) and (ii), and between case (i) and Langevin dynamics. 

PACS numbers: 05.40.Jc 



I. INTRODUCTION 

The static and dynamic behavior of polymers in the di- 
lute limit is well understood in the 2D and 3D limits, both 
of which have been studied theoretically with and with- 
out hydrodynamic interactions (HIs)^"'*. The theoretical 
predictions have been confirmed by experiment^^^. Ex- 
tensions have been developed to cover, for example, ef- 
fects of semidiluteness and chain stiffne ss, an d also these 
theories are well documented and testecP^'^'i^. 

Recent developments in promising micro and nanoflu- 
idic technologies^ have sparked a renewed intere st in the 
statics and dynamics of confined polymer^^EI!^ Classic 
theoretical works on the subject based on scaling theory 
esta blishe d the de Gennes and Odijk regimes of confine- 
mentliSLS. The former regime has a polymer confined at 
least in one dimension into a gap whose size L is smaller 
than the chain's radius of gyration in bulk, Rg, but much 
greater than its persistence length Ip. The de Gennes 
theory has been corroborated using the Edwards-Singh 
approach^*. Odijk developed a theory to describe the re- 
lationship between the planar size of the chain and 
L as L approaches Ip. Motivated by research on DNA, 
these regimes have been supplemented (with still exper- 
imentally unverified relations) to address more intricate 
regimes of confinement for different types of polymers^- 
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and followed by Monte Carlo (MC) simulations without 
chain-wall electrostatic interactionJsSEI]^ Most recently, 
Dai et al^^ have done extensive work without His on 
characterizing the transition to the de Gennes regime 
from weak conflnement and from the de Gennes to the 
Odijk regime. They found to follow a single power 
law with an exponent of about 1/4 as a function of 1 / L 
from the de Gennes regime up to the onset of the Odijk 
regim e in s litlike confinement. Two recent DNA experi- 
ment422i2aj have measured similar power laws and report 
exponent values of /? = 0.23±0.03 and 0.26±0.04, respec- 
tively. However, in another recent experimental study, 
Strychalski alW found iiy - {l/L)~°-^^ in nanofiuidic 
slitlike confinement close to and in the transition between 
the de Gennes and Odijk regimes. Moreover, other ex- 
perimental findings'^^ of polymer diffusivity in similar 
nanofiuidic slitlike confinement have contradicted predic- 
tions based on blob theory^^ and suggest the transition 
from the de Gennes to the Odijk regime to be broad. 

The disagreement and controversy between theory and 
experiments leaves room for numerical simulations that 
take into account the full gamut of His in such quasi- 
two-dimensional (q2D) systems where the polymer dy- 
namics is strongly restricted in one dimensiorP^. Particle 
dynamics is sensitive to hydrodynamic boundary condi- 
tions {e.g. porosity, slip or no slip on walls ), wh ich can 
dictate the physical behavior in the systenPSsH. Thus, 
the methodology must be chosen carefully for studies of 
confined suspensions in order not to neglect important 
features of stispensions dynamics due to limitations in 
the algorithnpS. Here, we take advantage of a recently 
developed hybrid Lattice-Boltzmann (LB) modepSl that 
provides a well-defined hydrodynamic radius and bound- 
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ary conditions for the solute particle^^ll 

In the present work, we investigate the transition from 
2D to 3D both for static and dynamic scahng of a sin- 
gle polymer confined between two rigid parallel plates. 
We consider two distinct types of confinement: (i) where 
both the solvent and the polymer are fully confined (q2D) 
and (ii) where the solvent is in 3D and only the polymer 
is constrained by a focusing potential. We quantify the 
static and dynamic scaling exponents as a function of the 
degree of confinement. Our method allows us to study 
the role of no-slip walls on the planar center-of-mass dif- 
fusion coefficient of the polymer. 

Unlike colloidal suspensions^, polymer dynamics in 
the q2D-transitional region with His has not been in- 
vestigated systematically as a function of the degree of 
confinement. It is not even known conclusively whether 
static equilibrium quantities are affected by His as they 
are coupled dynamically through fluctuations. Despite 
the pivotal role of His, polymer chains in q2D sy stems 
are often studied in equilibrium in the absence of HlJ^SEi] 
Full solutions to this complicated hydrodynamic problem 
are needed in order to provide understanding of funda- 
mental physics in microfluidic experiments^^ and future 
applications'^'^ . 

II. THEORY 

In the present work, we shall examine polymers con- 
fined between parallel plates a distance Lz apart, as 
shown schematically in Fig. [l] The confinement can re- 
sult in the polymers adopting an asymmetric configura- 
tion which we characterize using the tensor of gyration 

Qap, 

1 ^ 

QaP = - rcm,Q)(r„,/3 - Tcm^p), (1) 

n=l 

where TV is the degree of polymerization, a, /? £ {x, y, z} 
and the subscript cm refers to the center of mass of the 
polymer. The components of (Q) will be compared with 
the isotropic radius of gyration in the absence of confine- 
ment, Rg, which is related to the trace of Eq. Q and the 
sum of the tensor's eigenvalues by 

i?2 = (TrQ) = (A1 + A2 + A3) (2) 

in free space. The subscripts of the non-negative eigen- 
values indicate their respective magnitude, i.e. Ai > 
A2 > A3 > 0. The ensemble averages (Ai)/(A2) and 
(Ai)/(A3) characterize the anisotropy of the mean shape 
of the chain. To relate our simulations to experimen- 
tal imaging of the polymer by a microscope, we define 
the size of the polymer parallel, and normal, R±, to 
confining plates through 

3 3 
ii;|^(^A,(f-|e,.zp)); i?i^(^A,|e,.zp), (3) 

i=l i=l 



where are the three unit-normalized eigenvectors cor- 
responding to Ai and z is a unit vector perpendicular to 
the confining planes. We shall examine the polymer be- 
tween parallel plates separated by a distance Lz, which 
we use together with Eq. ([2| to define the degree of con- 
finement, C, as 

C=^, (4) 

where Rg is the free-space radius of gyration. As the 
Odijk regime with His is not feasible for our numerical 
method, we concentrate on the de Gennes scaling regime, 
I < C ^ Rg/lp, where Ip is the persistence length of 
the chain. For this range of C, we recapitulate how the 
polymer size parallel to the plates, Ru, should scale as a 
function of C. De Gennes proposecpl that between the 
plates, the self-avoiding polymer forms a set of sequen- 
tially connected blobs of linear size Lz each of which con- 
tains g monomers. The part of the contour contained 
within a blob behaves like a self-avoiding walk (SAW), 
i.e. Lz ^ crg'^ , where v = 0.588 is the asymptotic value 
of the 3D Flory exponent for an unconstrained SAW at 
infinite dilution"^^. Finally, the blobs behave like a 2D 
SAW for which reason the planar size of the confined 
chain scales as 

i?ll - Lz{N/gY^'' - LzN^^^ia/LzY'"'^'' 

Rg 

where we have used the Flory radius of gyration^^ Rg ^ 
aN'^ . The long-chain limit for the exponent {3 = vixyjv — 
1 w 0.276, where v^d = 3/4 is the asymptotic value of 
the Flory exponent in It is important to remember 
that in actual numerical simulations'^'^ with finite N the 
value of V (1^20) is larger than 0.588 (0.75) due to finite- 
size effects. 

The equilibrium size of a polymer is related to the 
static structure factor^, defined as 

1 ^ 

= ^ (exp[*k.(r™-r„)]). (6) 

m,n— 1 

The structure factor of an unconfined self-avoiding chain 
at infinite dilution scales as S{k) ~ k~^^^ in the scaling 
region kRg ^ 1 [k — |k|). These scaling results are 
subject to finite-size corrections from which asymptotic 
behavior can be extracted.'^ 

The dynamic version of Eq. (|6]) defines the intermedi- 
ate scattering function (or, the dynamic structure factor) 
5(k,i), 

1 ^ 

^^^^^'^^N ^ (exp[ik-(f„,(i + s)-f„(s))])^, (7) 

m,n— 1 
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where f (t) = r{t) — rcm{t)- The subtraction of the center- 
of-mass (CM) motion in Eq. ([t]) is important in order to 
extract scahng for the intramolecular dynamics only!^ 
Also, most close-range contributions to the sum in the 
analogue to Eq. Q in the laboratory frame are between 
consecutive monomers along the backbone, which results 
in scaling corresponding to semi-flexible chains in the lab- 
oratory frame of reference jiSiSMZl The scaling prediction 
for Eq. ^ i^ 

S{k,t) = S{k,0)F{k't), (8) 

where F{k^t) is a scaling function and z = 2+i'd/i' is the 
dynamic scaling exponent. Equation Q is expected to 
hold at infinite dilution for a long polymer at intermedi- 
ate lengthJS/c G (27r/i?g, 27r/(T), where cr is the size of the 
monomer. The quantity i^d is the scaling exponent asso- 
ciated with how the center of mass diffusion coefficient of 
the polymer, D, scales as a function of iV: Z3 ~ N~'^° . 
In 2D, it is known that — Q hi the presence of hy- 
drody nami c interactions as I? ~ log A'^ for which reason 
z = 2.121221 Jn 3D, y and vd are equal according to the 
Zimm model resulting in z = 3. Self-avoiding chains in 
the absence of hydrodynamic interactions have vd = ^ 
for dilute 3D systems and for all polymer concentrations 
in 2D leading to z « 3.7 (3D) and z = 10/3 (2D). 

We analyze the static and dynamic structure factors 
separately parallel to confining planes, S^^, and perpen- 
dicular to them, Sj_, which we define as 

5||(fc,t) = 5(k||,t) = S{k{e,+ey)/V2,ty, (9a) 
S±{k,t) ^ S{k^,t) ^ S{ke„t). (9b) 

We expect S'y to transition from the known 3D to 2D 
result with increasing level of confinement in the z direc- 
tion. The case of S± is not as clear at high degree of 
confinement and we do not expect it to exhibit scaling. 




FIG. 1. An illustration of our polymer model and coordinate 
system in the present work. The frame encloses a schematic of 
the structure of a 20-node composite monomer that provides 
hydrodynamic consistencj^^. In this work, the node count is 
30. The picture of the polymer was generated in vmeP. 



where the energy scale is set by e = A;b(300K), (T,nw = 
1.87 fj, fee is the Boltzmann constant and a — 1.5 nni. 
The interaction is cut off at 2^/^(Tmw ~ 2.1 cr. An anal- 
ogous potential is used at the top. We define the plate 
separation in Eq. Q as 

= (A^zAx - Cr,nw) - (Ax + Cr,nw)- (H) 



B. Polymer 

We use the finitely extensible nonlinear elastic (FENE) 
chairP^ to model a linear polymer as A^ monomers whose 
centers of mass (CMs) are connected by A^ — 1 bonds. To 
capture the excluded- volume effect, we apply the shifted 
and truncated 12-6 LJ potential between the CMs of all 
pairs of monomers. These two features can be written 
for a pair of particles i,j as 



III. MODEL 



A. Geometry 



The system is bounded by two parallel plates spanning 
a planar region [O,^^;] x [0,Ly] located at coordinates 
zq — Ax and N^Ax, where is a number of fluid lat- 
tice sites chosen to give a desired degree of conflnement. 
Figure [T] contains a schematic side view of the system. 
Periodic boundary conditions (PBCs) are applied in the 
X and y directions with L^. = Ly and they are chosen to 
be about lORg, which we find sufficient for eliminating 



finite-size effects in Section flV A 1 The polymer is con- 
fined between the plates by a shifted and truncated 12-6 
Lennard- Jones (LJ) potential normal to the planes. At 
the bottom wall this is 



Uir,,r,)^--kRl\og 



Rl 



+4e 



where rij = jr^ — rj|, the first term is the FENE bond 
present only for consecutive beads of the chain, Q{x) is 
the Heaviside step function and the second term is the 
repulsive LJ potential. The parameters of the model are 
the spring constant k = 32.6ecr~^ and the maximum ex- 
tension of a bond i?o = l.Scr. Each monomer i follows 
Newton's equation of motion 

mv, = - ^V,C/(r,,rj) +Fi(r,,Vi,u(r,,i)), (13) 



Z- Zq 



Z- Zq 



(10) 



where the mass of the unsolvated particle is m = 24 x 
10~24 g_ ^jj^g second term on the right-hand side is a local 
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force to describe frictional fluid-particle interaction, 



-ft d^x 



7ni(r)(v,+WjX(r- rj)-u(r))d^a;,(14) 

Bdt) ^ ' 

where is the monomer's angular velocity and V; its 
center-of-mass (CM) velocity. The local fluid velocity 
u(r) that enters Eq. ( 14 1 provides correct thermalization 
for the chain without having to add external noise to 
the particles when 7 is chosen appropriatel} P°P^l The 
composite monomer of volume Bi consists of 30 nodes 
arranged on a roughly spherical shell as shown in the in- 
set of Fig. [TJat a distance of i? = 0.7 nm about the CM 
coordinate r^. By modeling a monomer in this fashion, 
correct thermalization and consistency between different 
measures of the monomer's hydrodynamic radius is at- 
tained by choosing the coupling constant 7 as detailed in 
Ref.Eni 



C. Solvent 

1. Lattice-Boltzmann fluid 

Our solvent model reproduces equations for the mass 
and momentum conservation in a fLuid at the Navier- 
Stokes level that read 



dtp + da{pua) = 



(15) 



and 



dt{pUa) + df}{pUaUi3) ^ -daPaP + fa (16) 
+ dp (^{daUfi + dpUa - '^d^U^Sai3^ + (dyU^SalS^ , 

where p and Uq, are the fluid density and components 
of velocity, rj and C a-re the shear and bulk viscosities 
and Pa/3 is the fluid pressure. In this work, the pres- 
sure tensor is diagonal with linear dependence on den- 
sity, i.e. Pap = pVgSap + Sap, where Vs is the speed 
of sound and Sap is related to thermal fluctuation^. 
This can be viewed as an ideal gas equation of state 
or the first term in a Taylor expansion of the pressure 
about fixed density in which case is the isentropic 
compressibilitj0i]_ xhe force density exerted by the poly- 
mer appears through /„. Our standard LB fiuid algo- 
rithm reproduces Eqs. ( |15[) a nd ( 16 1 in the form typi- 
cal to most LB algorithmsP^ The shear viscosity in the 



model is 77 = p{t — At/2)i;^/3, where 



Ax /At is a 



lattice velocity, and ( = ri{5/3 — 3v'^/v'^)^l-n. this paper, 
T — 1.13 At and p = m/Ax^. The speed of sound Vg is 
chosen to be Vs — Vc/V^ (fs < fc is required for stabil- 
ity in LB algorithms). Thermal fluctuations in the fluid 
stress tensor Sap, and corresponding noise in higher mo- 
ments, maintain a constant temperature of T = 300 K.'^'^ 



Within the simulations themselves, computation is done 
in units of the lattice discretization Ax = 1.0 nm and 
At = 0.15 ps. The fluid boundaries are located half a 
lattice constant outside the MD walls, at z — Ax/2 and 
z — (Nz + 1/2) Ax, and they are implemented using the 
mid-grid bounceback rule''^ . is the number of lattice 
sites in the z direction. 



2. Langevin dynamics 

We also ran Langevin dynamics (LD) simulations of 
a chain of point particles in which the LB coupling was 
completely absent. We set the point monomer mass to 
2m. The fluid-particle interaction of Eq. ( 14) for the LD 



runs was simply replaced by the common combination of 
a frictional force and a Gaussian random force uncorre- 
lated in space and time. 



Fa = -(Va + P 

(P«(r„t)F«(r„t'))=2fcBTex 



R. 



(17a) 
(17b) 



,)S{t-t')6ap. (17c) 



The friction in the Langevin equation was set to ^ = 
67r?7_R, where R is the hydrodynamic radius of the 
monomers in the LB simulations. The time step in 
Langevin simul ation w as set to At = 35fs « 0.007tlj, 
where tlj = a^2m/(.. 



IV. RESULTS 
A. Statics 

The size and shape of a long polymer in a sli t has 
been studied recently using Monte Carlo simulationJ^SM]^ 
These studies had extensive sampling even for chains of 
several hundred monomers. However, hydrodynamic in- 
teractions were not included. We have simulated chains 
of = 32 to 96 monomers with (LB) and without (LD) 
hydrodynamics. Our study spans from weak confine- 
ment, C ~ 0.2, to the de Gennes regime, 1 < C ^ Pg/^v 
where 1^ is the persistence length of the chain. 



1. Monomer distribution 

We start by examining how monomers are distributed 
in the slit by looking at the in-plane distribution of 
monomers P\\ (r) about the chain's center of mass, where 

= (E„|rn,|| - r||,cml/^), and r„j| = (r„_^,r„_y) is 
the in-plane location of the nth monomer. The no-slip 
walls destroy momentum conservation in the z direction 
and reduce hydrodynamic finite-size effect from 1 /L 
to 1/L^, which is a consequence of t he fo rm of pairwise 
His (mobility tensor) in q2D system^'''^^ 
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We examine three cases. In all of them {N, C) = 
(96, 2.0) and the potential of Eq. ( [lO| confines the poly- 
mer between zq = Ax and N^Ax. In the first case, 
the LB fluid box size is {L^,Ly,L^) « (12, 12, 0.5)i?g 
(squares in Fig. [2]). No-slip walls for the fluid are in 
place at the top and bottom of the system. As poly- 
mers can be trapped experimentally in slit-like confine- 
ment either by walls acting on both the polymer and 
solvent or by a potential acting solely on the polymer 
{e.g. via an optical trap), we also study the second 
type of confinement. Therefore, in the second case, 
the LB fluid box is of size {12Rg)^ (circles) and has 
periodic boundary conditions (PBCs) in all directions. 
This changes the form of His from q2D (squares) to 3D 
(circles). We also include Langevin dynamics (LD) in 
{Lx,Ly,Lz) « (12, 12, 0.5)i?g as the third case (trian- 
gles). The resulting (r) for each case is shown in Fig. [21 

The planar radius of gyration, , is related to (r) 
via 



This formula gives Rii 



1.19R 



\\,bulk 



(circles) and 1.41i? 



r^P\\ (r) dr. (18) 

= 1.18i?|j,fc„;fe (squares), 
(triangles) for the 



\\,bulk 



three cases mentioned above. These values, and the 
plot of P\\{r) shown in Fig. [2j indicate the monomer 
distribution for the 3D-LB fluid to be much like that 
for the q2D-LB fluid (squares vs. circles). The LD 
chain extends 23% more (i?g(LD)/i?g(LB) sa 1.23) at 
this level of confinement. Albeit weak, we note that the 
joint effect of no-slip walls (see text and Fig. [3] below) 
and the form of His in q2D is observable for a poly- 
mer at a level of confinement C ~ 2 as a slightly fat- 
ter tail (squares vs. circles in Fig. [2|. Even though the 
data is smooth, the slightly fatter tail in q2D than in 
3D cannot be attributed directly to an effec tive pair- 
wise antidrag in q2D between monomer^^MIl since the 
slight difference occurs around r w 1.5i?||_buik ~ 2.3Lz. 
This is not in the limit r ^ where q2D interactions 
should hold for blobs in de Gennes' theory. Going from 
an unconfined fluid (3D-LB, circles) to q2D (squares) 
changes the form of long-range interparticle interaction 
from l/r to l/r^, which also reduces finite-size effects 
due to self-interactions across periodic boundaries. We 
found finite-size effects to be negligible by simulating 
{N,C) = (32,1.6) for = Ly 
saw no discernible changes in Py (r) between them 

We will also later compare our results to LD simu- 
lations. A precise quantitative comparison between LD 
and LB simulations can be problematic even if the MD 
parameters are equal as the presence of His and no-slip 
boundary condition modifies the level of confinement in 
the z direction. This can be seen in the probability dis- 
tribution P±{z) (normafized as J^"'^^ P±{z) dz — 1) of 
monomer position between the plates shown in Fig. [3] 
Monomers on the Langevin chain (triangles) are able to 
get closer to the MD wall than those in a q2D-LB (cir- 
cles) or a 3D-LB fluid (squares). The no-slip walls in q2D 



lOPg to 20Pg and 



: (12,12,0.5)Rg, q2D-LB 
: (12Rg)^ 3D-LB fluid 
A : (12,12,0.5)Rg,LD 




,bulk 



FIG. 2. Probability distribution Py (r) of the planar monomer- 
CM distance r = {J2„ lr„,|| -r||,cm|/A^) for (iV, C) = (96,2.0). 
The abscissa in the figure is scaled by the planar free-space 
radius of gyration of the chain. The legend refers to the LB 
and LD box sizes, see text for details. 
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FIG. 3. Half of the symmetric probabiHty distribution P±{z) 
of monomer positions between confining plates as a function 
of scaled distance from the center of the channel for a 96-bead 
chain for q2D-LB (squares) , 3D-LB (circles) and Langevin dy- 
namics (triangles) at C = 2.0. The dashed line is located at 
the cutoff (2^''^(Tniw from the wall) of monomer-wall interac- 



tion and the dotted line is at (Tmw from the wall, see Eq. (|10| . 
The thinner tail of the LD distribution as compared to q2D- 
LB tells how for a given C, the two types of dynamics yield 
effectively different degrees of confinement. 



provide additional repulsion compared to the 3D-LB fluid 
(circles vs. squares). 



2. Radius of gyration 

Figurefflshows the normalized planar, and perpen- 
dicular, /Tl, radius of gyration as a function of 1/C for 
both q2D-LB (solid symbols) and LD (hollow symbols) 
simulations. We flnd signiflcant deviations in R± and P|| 
from bulk behavior for 1/C < 2. Linear flts to the loga- 
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FIG. 4. Normalized (a) planar and (b) perpendicular size of 
the polymer in log-log scale as a function of the inverse of the 
degree of confinement for LD simulations without His (hollow 
symbols) and q2D-LB simulations with His (solid symbols). 
The MD confining plates according to Eq. ( |10[ | are identical in 
all cases. Data points are for A'^ = 32 (large squares), 64 (large 
circles) and 96 (up triangles) so 1/C = const, may correspond 
to different Rg-L^ combinations. The symbols corresponding 
to 1/C > 2 are excluded from the fits whose results are shown 
as straight lines for q2D-LB data and dashed lines for LD 
data. 



rithm of the data yield scaling exponents /3 = 0.24 ± 0.03 
for q2D-LB and 0.21 ± 0.03 for LD. The scaling of the 
perpendicular component R± in Fig. [4] is close to linear 
indicating that for 1/C = L^/Rg < 2 the chain takes up 
all the available space in the z direction independent of 
the type of dynamics. The LD chain has a larger rel- 
ative extension in the parallel component compared 
to the LB chain. We attribute this to the absence of 
solvent-mediated correlations and the missing resistance 
to monomer movement perpendicular to the backbone 
of the chain. Both data sets (q2D-LB and LD) contain 
points for iV = 32, 64 and 96. In the perpendicular direc- 
tion there is very little difference between q2D-LB and 
LD results as the chain takes up all available space in 
that dimension. However, the approach to unity is much 
slower than in the parallel direction. 

Quantitative comparisons to the theoretical predic- 
tion of Eq. ([5]) have been performed in recent years 
both experimentally^^"^^'^®' and by computer simulations 
witlF^ and without HIs^^'SJI ^ote that some pub- 



lications define C in Eq. ^ as i?||,buik/-^z instead of ac- 
cording to Eq. This ooes not affect the scaling ex- 



ponent as the relation i?||^buik — 
space due to isotropy. 



holds in free 



We comment on the results of computer simulations 
first. The value of /3 depends on the 3D Flory expo- 
nent, which itself depends on A'' due to finite-size ef- 
fects for the range of N considered here. We have ob- 
servecPSl values from i^{N = 32) = 0.628 ± 0.004 to 
,y{N = 96) = 0.595 ± 0.003 in bulk with His, from which 
we could extract the value i' = 0.586 ± 0.005 for the 
asymptotic iV — >■ 00 limit. Dai et al. used N — 101 to 
401-bead chains at C £ (0.1,37) in their Monte Carlo 
studj-^of slitlike confinement. In order to define bounds 
for the de Gennes regime, they chose € [2/p, i?j|.buik]j 
which kept deviations between data points and the line of 
best fit below 6% for C £ (0.6, 3.1). These criteria yielded 
power law exponents of I3{N = 101) = 0.236 ± 0.023 to 
I3{N = 401) = 0.255 ± 0.003, which translate to v{N ^ 
101) = 0.607 ± 0.011 and i/(7V = 401) = 0.598 ± 0.001 
with 1^20 = 3/4 using Eq. We note that also the value 
of i^2D is larger than 3/4 for finite N, which, if used in 
calculating iy{N), would only increase the value of i^{N) 
based on a given /?. Dai and coworker^^Sl -went further 
by studying the shift from the de Gennes to the Odijk 
regime, but this was carried out using the effective width 
of the chain as the scaling parameter which would corre- 
spond to C > Rg/{2lp) in our model. However, the width 
parameter is not included in the our model for which 
reason we cannot access the Odijk regime in the present 
work. Micheletti and Orlandini^i found /3 = 0.27 ± 0.05 
for C e [2,5] and N e [360,480] in their MC work on 
semi-fiexible chains. Our LD simulations for 32, 64 and 
96-bead chains yield /3 — 0.21 ± 0.03 (hollow symbols in 
Fig.|4]). Also, our resuh /3 = 0.24±0.03 based on q2D-LB 
data agrees well with the result /3 « 0.24 from an early 
Brownian Dynamics study^^. Guided by these numbers, 
it is fair to state that computer simulations of coarse- 
grained flexible and semi-flexible polymers both with and 
without His, considering the finite value of N, support de 
Gennes' theory of Eq. ([5]) if 1 < C <C Rg/lp- However, 
due to the limited range of values of C that span less 
than a decade in the present and previous simulations, 
it is impossible to conclusively verify the existence of the 
de Gennes scaling regime. 



Two recent DNA experiment j^^ l ^'^l report values of 
/3 = 0.23 ± 0.03 (2.25 < C < 7.8) and 0.26 ± 0.04 
(0.4 < C < 1.7 ), respectively. However, other recent ex- 
perimenti^mil probe levels of confinement up to C w 20 
and C ~ 132 and suggest lower values of /3 = 0.17 (based 
on figure 4a in Ref. [26|) and /? = 0.16 ±0.01, respectively. 
We note that the onset of the theoretical Odijk regime 
[Lz ~ Ip, where Ip is the chain's persistence length) cor- 
responds to 12 < C < 20 for double-stranded A-DNAI^I!. 
Such high values of confinement are beyond the scope of 
the type of model we consider here which loses its phys- 
ical relevance in the Odijk regime. 
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FIG. 5. Both the (a) planar and (b) the perpendicular compo- 
nent of the static structure factor show differences when mea- 
sured in q2D-LB (solid lines) and 3D-LB simulations (dashed 
lines). The polymer conta ins 32 beads. The insets show the 
differences defined in Eq. ( 19 1 at different C as a function of 
ka. The lines in the insets correspond to C = 1.6 (solid), 0.8 
(dashed), 0.4 (dotted) and 0.3 (dot-dashed). 



Static structure factor 



FIG. 6. The normalized difference in the (a) planar and (b) 
perpendicular static structure factor between q2D-LB and 
LD simulations at different degrees of confinement (Eq. Q). 
Structural differences can be seen at values of ku correspond- 
ing to the monomer size {ka ~ 6.3), the maximum bond 
length ka — 4.2, the spacing between plates a,t C — 1.6 
(solid line at ka ~ 3.5), the distance between next-nearest 
neighbor monomers {ka « 2) and twice the radius of gyration 
{ka = 0.7, ...,0.9). 



Y. von Hansen et al. found that as a polymer's mean 
distance to a wall is decreased, there is a slow transition 
from Zimm to Rouse dynamics in terms of hydrodynamic 
screening efFectJ^. Their observation leaves room for the 
possibility that even static quantities such as the static 
structure factor could be affected by the presence of con- 
fining walls. We investigate this in Fig. [Sj which shows 
the planar {S\\{k), Eq. (j9^)) and perpendicular {Si_{k) 
Eq. ^p)) components at different levels of confinement 
for q2D-LB fluid confined between two no-slip walls (solid 
lines in Fig.jsJ^a) and (c)) and 3D-LB fluid (dashed lines 
in Fig. [5ja) and (c)). Two parallel plates confine the 



chain in both cases according to Eq. (10 1. A comparison 



between the q2D-LB and the 3D-LB fiuid (a box of size 
(52Ax)'^ with PBCs) allows us to ascertain whether an 
effect due to direct momentum absorption by the no-slip 
walls in the z direction can be seen in the static structure 
factor. Figures [sjja) and (b) show the differences 



As(fc) 



S'q2D,LB(fc) - 5'3D,LB(fc) 



3D, LB 



(fc) 



(19) 



at different levels of confinement in the plane of the walls 
(II) and perpendicular to it (_L). Differences up to 4% (||) 



to 24% (_L) are visible at strong confinement (C = 1.6, 
solid lines), and at the opposite end (C = 0.3, dot-dashed 
line) , the maximum differences in the planar and perpen- 
dicular data are 4% and 7%. These observations sup- 
port the finding that i?g,_L approaches the unconfined 
limit more slowly than i?g y in Fig. 4j The no-slip wall 
has an observable effect on chain statistics in the per- 
pendicular direction at C = 0.8 and C = 1.6 for for 
ka > 1.5. The largest deviations in the parallel direction 
occur at C = 1.6 for 0.25 < fccr < 1.0, which corresponds 
to lengths {1.7Rg,7Rg). However, these are not statis- 
tically significant as the maximum statistical errors for 
that range of ka in the differences are 0.06 at C — 1.6 
(solid line) and 0.04 at C = 0.8 (dashed line) in the inset 
of Fig.^a). 

Next, we shall ascertain whether the static structure 
factor shows any differences due to the presence of hydro- 
dynamic interactions. This is done by comparing results 
from q2D-LB and LD simulations, which we comment on 
first. The monomer structure in LB is stiffer due to the 
composite nature of monomers of Fig. [l] 

Figure[6]shows differences in the static structure factor, 
but between q2D-LB and LD simulations, i.e. As(fc) = 



8 



('S'q2D,LB(fc) — S-^uik))/ S'LD(fc)- The planar data in panel 
(a) tell that the Langevin chain frequents configurations 
for which 1.3 < ka < 3 more than the LB chain (max- 
imum statistical error in As||(fc) is 0.04 for C — 1.6 
and 0.03 for C — 0.8) and the situation is reversed for 
0.4 < fccr < 0.9 (maximum error 0.09 for C = 1.6 and 0.05 
for C = 0.8). The former range of ka corresponds to Rg 
and the latter to 2Rg. The differences become smaller as 
C tends closer to zero (dotted and dot-dashed lines) as is 
expected for chains in bulk. The perpendicular data at 
C = 1.6 in panel (b) reveal the local stiffness of bonds in 
the LB fluid (fccr = (27r/(2.25 nm))cr « 4.2). Again, as C 
tends to zero, the difference in s±{k) becomes smaller. 

With the possibility of scaling in the form of a power 
law, we plot —l/i/ = d[log S (k)] / d[\og k] as a function of 
k in Fig. [7) We entertain the questions if and which scal- 
ing exponent the planar static structure factors exhibit. 
In approaching the 2D limit, C would tend to infinity and 
1/iy = 4/3 and in the 3D limit, C ^ and 1/;^ w 1.7 for 
a long chain in the scaling region. Figure [7] has two pan- 
els showing —1/iy for N — 96 extracted from the planar 
structure factor, S'||(fc), based on q2D-LB and LD simu- 
lations. The LB data have a decay exponent larger than 
1/v = 1/0.6 in the anticipated scaling region whereas 
Langevin simulations in Fig.[7](b) follow the line —1/0.75 
corresponding to 2D-like scaling for 0.8 < fccr < 2 es- 
pecially at the highest degree of confinement (C = 2.0, 
hollow squares in (b)). That is, the static structure in 
the presence of His suggests 3D dynamics whereas strong 
confinement in LD produces effectively 2D dynamics for 
the same value of C . The q2D-LB data comes closer to 
the —l/v2D line for C = 3.1 (not shown), which is the 
most confined case we can access. The monomer size 
corresponds to fccr = 2it and the radius of gyration to 
fccr = 0.8. Moreover, the unconfined limit C — is ap- 
proached differently by the two types of dynamics: LB 
data do not exhibit clear 2D scaling in the anticipated 
scaling region, but indicate 1 /v to decrease from about 2 
to 1/0.6 as C — >■ whereas LD data transition from 2D 
{l/v = 1/0.75) to 3D scaling (1/0.6) as C ^ 0. 



B. Dynamic scaling 

Equilibrium polymer dyna mics in bulk is described by 
the theory of dynamic scaling!^. The scaling prediction 
of Eq. (Is]) in the presence of His can be tested by plotting 
S{k,t)/S{k,0) versus (fccr)^(</TLj) for different fc, which 
should all collapse onto a single curve, i.e. F of Eq. ([s]), 
with 2: = 3 in three dimensions and z = 2 in two di- 
mension^2ll xhe corresponding theoretical values in the 
absence of His are z — 3.7 and z — 3.3§^. 

Since momentum conservation is broken in the z di- 
rection in our q2D system, one cannot expect either re- 
sult to hold as such. It is therefore interesting to see 
what the transition between 3D and 2D scaling behav- 
ior looks like. As in the static case, we treat the xy 
plane and the perpendicular z direction separately. We 
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FIG. 7. Exponent —l/f based on the relation S'||(fc) ~ k~^'''^ 
for a chain at levels of confinement C — 2.0 (squares), 0.9 
(circles), 0.5 (up triangles), 0.3 (down triangles) and 0.2 (di- 
amonds). Panels are for a 96-bead chain from (a) LB and 
(b) LD simulations. The top horizontal line indicates the 2D 
scaling exponent —1/1^20 = —4/3 and the bottom line the 3D 
exponent — I/i^sd = —1.7. 



find 5|| (fc, t) for = 32 with His to exhibit values of z 
that are close to the known 2D {z — 2) result in strong 
confinement and approach the 3D result [z — 3) as the 
level of confinement decreases. Evidence of this transi- 
tion is provided in Fig. |8j which shows scaling collapses 
for ^[[(fc,^) in which z increases from 2.2 at C = 1.6 
to 2.7 at C = 0.3. The C ^ limit was studied in 
Ref.inniand gives z = 3. The collapses are plotted as func- 
tions of (k^ty^^^ as it has been shown theoretically that 
log i^(a:) x^/^ for k'^k^Tt/ {<ot:vi) > 1 in free spacfPSl. 
Since the right-hand side of the inequality is at least 10 
for our parameters, we opted to use this form in the scal- 
ing plots. When probed perpendicular to the wall using 
k = fce^, 5'_L(fc, i)/S'j_(fc, 0) do not collapse for any z at 
C — 1.6 i.e. at plate separation Rgl'2- (not shown), but 
at C = 0.3 (plate separation of ~ 3Rg) plotted in Fig. |9j 
we do observe a reasonable data collapse for A^ = 32 with 
a scaling exponent of z = 3 indicating 3D behavior. As 
Langevin simulations for long chains are feasible, we pro- 
duced scaling collapses of (fc, t) for A = 96 at C = 1.8. 
We found the chain to scale very clearly with the 2D scal- 
ing exponent of z = 2 + 1/0.75 3.33 for fccr € (1.4, 2.8) 
(not shown). 
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FIG. 8. Dynamic scaling for A'^ = 32 in the presence of His along k — k{ex + By)/\/2, ka € (1.8,2.7), i.e. parallel to the 
confining walls. As the level of confinement decreases from (a) C = 1.6 to (b) C — 0.4 and to (c) C = 0.3, the scaling exponent 
increases from z = 2.2 ± 0.1 to z — 2.7 ± 0.1. 
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FIG. 9. For plates spaced C — 0.3 apart, the dynamic struc- 
ture factor in LB fluid for A'^ = 32 exhibits 3D scaling behavior 
in the direction perpendicular to the plate. 



C. Diffusion 

Recent experiments by Strychalski et al-^^ and by 
Tang et al^ have shown the degree of confinement to 
be reflected in the chain's planar diffusion coefficient as 
~ (0.4 < C < 25) and (0.3 < C < 21) in 

aqueous solution. Moreover, their respective findings 
a sa 0.48 ± 0.05 (as averaged over Strychalski's results 
for different chain lengths) and a = 0.49 zfc 0.01 are in 
agreement with previous experimentd^^'^, but in con- 
tradiction with the predictions of blob theory that gives 
a = However, Usta et alF^ iound agreement with 



blob theory in slit-like confinement using their Lattice- 
Boltzmann scheme. Their LB scheme differs from ours 
in that we need not add external Langevin noise to our 
chain as we have consistently calibrated the hydrody- 
namic radius of our monomers and the coupfing between 
the monomers and the LB fiuicP^. 

We define the in-plane center-of-mass diffusion coeffi- 
cient of the chain through the mean squared displacement 



as 



Dn = 



\t- 



lim 

-to|-i-oo 



((r||,cm(i) 



(20) 



where r||,cm is the 2D vector (xcm, J/cm)- Figure 10 shows 
D|| in log-log scale as a function of 1/C for Jy— 32, 
64 and 96. In general, we observe a reduction in the 
diffusivity as the level of confinement is increased and 
for a given value of C, the longer a chain is the slower it 
diffuses. 

We find a single power law exponent to fit our log 
vs. log (1/C) data nicely over a decade starting from 
1/C = 0.5. However, by including data points in the 
region 1/C < 0.5, we get D„ - C-°-i&±o.Q4 ^j. = 32 

(solid squares), ~ c-°-^^^°-"'^ for = 64 (sohd triangles) 
and ~ C7-o.60±o.o4 iv = 96 (solid circles). Linear 
interpolation suggests (7-"-66±o.05 j^j. jy _^ which 
agrees with de Gennes' theory. 

The parallel-plate arrangement is suitable for testing 
the effect of direct momentum absorption by the no- 
slip walls in the z direction on a dynamic quantity such 
as the chain diffusivity. Interestingly, once the no-slip 
walls are removed (MD walls are still kept in place) and 
the LB fluid is then contained in a periodic box of size 
{52Ax)^, for = 32 (solid diamonds) becomes in- 
dependent of C. Even at 1/C « 3.6, the effect of q2D 
His and the non-conserved z momentum on is still 
14% (£'||.q2D/i'||.3D « 0.86). The LD simulations for 
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FIG. 10. The scaled planar center-of-mass diffusion coefficient 
as a function of 1/C in q2D-LB for A'' = 32 (solid squares), 
N = 64: (solid triangles) and A'^ = 96 (solid circles), in 3D- 
LB for = 32 (solid diamonds) and LD for iV = 96 (hol- 
low circles). The ffts to q2D-LB data give ~ c-°-'^^^°-°^ 



(sohd line, N = 32) 

-0.60±0.04 



c- 



(dashed line, = 64) 
and ~ C'-" '="=^" "'* (dotted line, N = 96). The errors are of 
the size of the symbols. The dot-dashed line with a slope of 
2/3 corresponds to the prediction based on blob theory. The 
scaling factor Do = k-BT{N(^)~^ is the diffusion coefficient for 
A'^ = 96 based on Rouse dynamics. 



= 32 (not shown for clarity) and 96 (hollow circles) 
are also unaffected by a change in C, which together 
with the LB simulations reveals the importance of the 
no-slip boundary condition with respect to the experi- 
mental observations. We note that even at the strongest 
confinement (1/C « 0.3), £)||(LB,iV = 96)/L»||(LD, TV = 
96) w 1.3 and thus, we concur with the conclusions by 
von Hansen et al^ regarding the broadness of the tran- 
sition from Zimm to Rouse dynamics in the vicinity of a 
wall. 



V. CONCLUSIONS 

We have simulated solvated coarse-grained polymers 
between two parallel plates of varying separation with 
(The Lattice-Boltzmann Method) and without (Langevin 
Dynamics) hydrodynamic interactions. As our numerical 
method is not applicable for studying the Odijk regime, 
we have only used plate separations that correspond to 
a regime where the de Gennes blob theory should hold. 
The parallel-plate system is unique due to the form of 
pairwise hydrodynamic interaction^^!! ^nd as the poly- 
mer moves almost in two dimensions when the plate sep- 
aration is comparable to or smaller than the polymer's 
free-space radius of gyration. For fixed MD confinement, 
we found hydrodynamic interactions to change the pla- 
nar monomer distribution about the chain CM slightly 
as the dimensionality of the fiuid was changed from q2D 
to full 3D. For levels of confinement C = 0.3 — 2.0, the 
Flory exponent based on the static structure factor at- 



tained values that varied between the free-space 2D and 
3D values depending on the length scale in question. 
The planar diffusion coefficient for the CM of the chain, 
D|| , decreased as a function of increasing confinement 
in agreement with experiments^^ ^® ''^ . For the lim- 
ited range of C accessible in our simulations, this sug- 
gests the scaling behavior to be due to the presence of 
no-slip walls alone and electrostatic interactions are not 
needed to account for the observed reduction. However, 
the value o f the scali ng exponent seemed to agree with 
experiment^^milllllll for short chains and asymptotically 
with de Gennes' prediction. As experiments are carried 
out close to the long-chain limit and are unaffected by 
finite-size effects, more studies are needed to ascertain 
the role of electrostatics in the scaling of the planar diffu- 
sion coefficient. Langevin dynamics (LD) simulations do 
not capture the reduction in Z)|| as a function C for which 
reason they are not a suitable tool for the study of dy- 
namic phenomena in a confined environment. However, 
the LD results are qualitatively similar to confinement 
imposed only on the polymer, but not the solvent. We 
have also characterized the differences in the static struc- 
ture factor between different types of confinement and 
dynamics. Finally, we observed a continuous transition 
in the polymer's dynamic scaling exponent from its value 
in 2D to the one in 3D as a function of decreasing level 
of confinement in the de Gennes and weak-confinement 
regime. 
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